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Abstract 

We evaluate the pseudoscalar-meson octet-baryon coupling constants and the corresponding ax- 
ial charges in eight channels (7rAW, 7t£E, 7tAE, KAN, KI.N, nEE, KAE and KEE) in lattice 
QCD with two flavors of dynamical quarks. The parameter a = F/(F + D) representing the 
SU(3)-flavor symmetry is computed at each u,d-quark hopping parameter and at the flavor-SU(3) 
symmetric point where the three quark flavors are degenerate at the physical s-quark mass. In 
particular, we obtain a - 0.395(6) at the SU(3) symmetric limit. The quark-mass dependences 
of the coupling constants are obtained by changing the u- and the li-quark masses and we find 
that the SU(3)-flavor symmetry is broken by only a few percent at each quark-mass we employ. 
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1. Introduction 

It is long after the Quantum ChromoDynamics (QCD) was established as the underlying 
theory of the strong interaction in hadrons or nuclei. On the other hand, QCD-based description 
of hadron physics has not yet been successful due to the strong coupling nature of QCD. The 
low-energy QCD dynamics is far from the perturbative regime, and nonperturbative analyses are 
essentially needed. Among potentially workable methods, lattice QCD has been established as 
a powerful first-principle method for QCD. Lattice QCD has been well developed and almost at 
the mature level. Nowadays, the physical point can be directly accessed with lattice QCD, and it 
can cast light on hadrons or hadronic interactions directly from the viewpoint of QCD. 

In hadron physics, meson-baryon coupling constants are one of the most important ingredi- 
ents. Baryon are the building blocks of our world, and meson-baryon coupling constants provide 
a measure of baryon-baryon interactions in terms of One Boson Exchange (OBE) models. In phe- 
nomenological potential models, the meson-baryon coupling constants are determined so as to 
reproduce the nucleon-nucleon, hyperon-nucleon and the hyperon-hyperon interactions in terms 
of, e.g., OBE models. It is then an important issue to determine such coupling constants directly 
from QCD, the underlying theory of the strong interactions. Lattice QCD is the only method 
that provides a first-principle calculation which can be directly compared with experiments, and 
it serves as a valuable tool to determine the hadron couplings in a model-independent way. 

Meson-baryon coupling constants directly reflect the flavor structure in hadron sectors. In the 
SU(3)-flavor [SU(3)f] symmetric limit, one can classify the pseudoscalar-meson-octet-baryon 
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coupling constants in terms of two (undetermined) parameters: the ttNN coupling constant and 
the a = F/(F + D) ratio H; 



gnNN = g, girJX = Iga, g nA z = — -g(l - O), gKAN ~ —g(l+2a), 

V3 V3 
gKW = g(l - 2a), gxss = -g(l - 2a), gKAB = ~p.g(^a — 1), gKiz = -g. (1) 

This systematic classification governs all the meson-baryon couplings in the case of SU(3)f 
symmetric limit. However as we move from the symmetric case to the realistic one, SU(3)f 
breaking inevitably occurs as a result of the mass difference between u, d- and s-quarks. The 
broken symmetry no longer guarantees a pattern for the meson-baryon coupling constants, and 
the breaking pattern can be determined only from the underlying QCD dynamics. They therefore 
should be individually calculated directly based on QCD. In this work we extract the coupling 
constants g„ m , g nIS , g nAZ , g KAN , g K1N , g nSS , g Khz and g K13 (denoted by g MBB > hereafter) by 
means of lattice QCD with two flavors of dynamical quarks. 



2. Meson-Baryon coupling constants and axial charges 

The pseudoscalar current matrix element is written as 

(!B(p)\P(0)\!B'(p')) = g P (q 2 )u(p)iy 5U (p'), (2) 

where gp(q 2 ) is the pseudoscalar form factor, — p' — p^ is the transferred four-momentum 
and P(x) = ^r(x)iysj-t//(x) is the pseudoscalar current. We compute this matrix element using the 
standard techniques on the lattice 10. S B S S 0> S • The baryon interpolating fields employed 
in this study are given as 

Tj N ( x ) = f abc [u ra (x)Cy 5 d b (x)}u c {x), rfr(x) = ^ bc [s Ta (x)CysU b (x)]u c (x), (3) 

77a (x) = \[u Ta (x)Cy 5 s b (x)]d c (x) - [d Ta (x)Cy 5 s b (x)]u c (x) + 2[u Ta (x)Cy 5 d b (x)]s c (x)} , (4) 

Vo 

77=(x) = e abc [s Ta {x)Cy 5 d b {x)V{x\ (5) 

where C = y\yi and a, b, c are the color indices. Calculating gp(q 2 ) (gp on the lattice) at the 
momentum transfers q 2 a 2 = n(2njL) 2 (for the lowest nine n points), where L is the spatial extent 
of the lattice, we obtain the meson-baryon form factor gMBB'iq 2 ) via the relation, 

L f 2n G M gMBB'iq 2 ) ft ~. 

g P (q ) = — 5 — r~ ' (6) 
m M - r 

assuming that the pseudoscalar form factors are dominated by the pseudoscalar-meson poles. 
Here Gm = (vac|P(0)|M) is extracted from the two-point mesonic correlator (P(x)P(Q)). Finally 
we extract the meson-baryon coupling constants gMBB 1 — gMBB'(0) by means of a monopole form 
factor: 

A 2 

gMBB'iq ) = gMBB'— 2 (') 

MBB' 1 

We employ a 16 3 x 32 lattice with two flavors of dynamical quarks and use the gauge config- 
urations generated by the CP-PACS collaboration [9] with the renormalization group improved 
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gauge action and the mean-field improved clover quark action. We use the gauge configura- 
tions at /? = 1.95 with the clover coefficient csw = 1-530, which give a lattice spacing of 
a - 0.1555(17) fm (cr 1 = 1.267 GeV), which is determined from the p-meson mass. The 
simulations are carried out with four different hopping parameters for the sea and the u,d valence 
quarks, K sea ,^, = 0.1375, 0.1390, 0.1400 and 0.1410, which correspond to quark masses of ~ 
150, 100, 65, and 35 MeV, and we use 490, 680, 680 and 490 such gauge configurations, respec- 
tively. The hopping parameter for the s valence quark is fixed to r 5 . = 0.1393 so that the Kaon 
mass is reproduced 19J], which corresponds to a quark mass of ~ 90 MeV. We employ smeared 
source and smeared sink, which are separated by 8 lattice units in the temporal direction. Source 
and sink operators are smeared in a gauge-invariant manner with the root mean square radius of 
0.6 fm. All the statistical errors are estimated via the jackknife analysis. 
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Figure 1 : LEFT: a a = F/F + D ratio extracted by a global tit of eight meson-baryon coupling constants is plotted as a 
function of mj. RIGHT: The SU(3)-breaking parameter [ 7, 8] extracted by a global fit of eight meson-baryon couplings. 
The empty circles in both figures denote the data at the SU(3)f symmetric limit. 

In the SU(3)-flavor symmetric case, where k"'^ = k s , = 0.1393 and the SU(3) relations are 
exact, all the couplings are well reproduced with a = 0.395(6), which is obtained by a global 
fit of the lattice data using Eq.(Q~|i. The obtained value of a is consistent with that in the SU(6) 
spin-flavor symmetry (a — 2/5) ifioll . which is the symmetry based on the nonrelativistic quark 
model. As we move to the chiral limit, the mass difference between u, d- and s-quarks gets large, 
and the fitted a gradually decreases toward the chiral limit. The quark-mass dependence of a can 
be found in the left panel in Fig. [TJ We here define and evaluate the SU(3) breaking parameter 
<5su(3), which represents to what extent the SU(3) relations are broken. (See Refs.fl7|, ISD for the 
detail.) The breaking parameter <5su(3) is plotted in the right panel in Fig.Q] Interestingly enough, 
<5su(3) remains 0(1)% in the quark-mass range we consider, which suggests for the pseudoscalar- 
meson couplings of the octet baryons that SU(3)f is a good symmetry. We try a quadratic fit of 
<5su(3) and find that 5su(3) = 0.072(16) in the chiral limit. The SU(3)f breaking effect seems to 
appear in a rather than in SU(3)f relations. 

We also extract the corresponding eight axial charges (g A w< Sajs, 8a,a^ 8a,an, gA,ZN, 8a,ss, 
£a,ae and gAjs)> an d repeat our analysis for them. Fig.|2]shows a and <5su(3) extracted by global 
fits of the eight axial charges. The SU(3) parameter a is found to be 0.390(2) in the SU(3) limit, 
and we observe essentially the same tendency as meson-baryon couplings. 



3. Summary 



We have evaluated the pseudoscalar-meson-octet-baryon coupling constants as well as the 
corresponding axial charges in eight channels (nNN, nYL, nAL, KAN, KI.N, 7rSH, KAE. and 



3 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

mk [lattice unit] mil [lattice unit] 



Figure 2: LEFT: a a = F/F + D ratio extracted by a global fit of eight axial charges is plotted as a function of mj. 
RIGHT: The SU(3)-breaking parameter extracted by a global fit of eight axial charges. The empty circles in both 
figures denote the data at the SU(3)f symmetric limit. 

KEE), in two-flavor lattice QCD with the hopping parameters K sea , k"^ = 0.1375, 0.1 390, 0. 1400 
and 0.1410, which correspond to quark masses of ~ 150, 100, 65, and 35 MeV. Our results 
indicate, for both of pseudoscalar-meson couplings and axial charges of the octet baryons, 

1. SU(3)f parameter a = F/(F+D) is 0.395(6) for the meson-baryon couplings (0.390(2) for 
the axial charges) in the SU(3)f symmetric limit, which is consistent with the prediction 
from SU(6) spin-flavor symmetry (a = 2/5). 

2. SU(3)f relations hold well, which is broken by only a few percent (at least) in the 35 MeV 
to 150 MeV range of the light quark masses. SU(3)f is hence a good symmetry in the 
quark-mass range we consider. 

3. The parameter a gradually decreases towards the chiral limit, and SU(3)f-breaking effect 
appear in a rather than SU(3)f relations themselves. 
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